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' Abstract 

In this paper the arising of Gribov copies both in Landau and Coulomb gauges in regions with 
cm : non-trivial topologies but flat metric, (such as closed tubes S 1 x D 2 , or RxT 2 ) will be analyzed. 

Using a novel generalization of the hedgehog ansatz beyond spherical symmetry, analytic examples 



\& ' of Gribov copies of the vacuum will be constructed. Using such ansatz, we will also construct the 

elliptic Gribov pendulum. The requirement of absence of Gribov copies of the vacuum satisfying the 
strong boundary conditions implies geometrical constraints on the shapes and sizes of the regions 
. with non-trivial topologies. 

m 

1 Introduction 

X 

One of the most important sectors of the standard model is Yang-Mills theory which describes 
QCD and the electro-weak theory. The degrees of freedom of a gauge theory are encoded in the 
connection [A^) 11 , which is a Lie algebra valued one form. The action functional is invariant under 
gauge transformations 

A„ -> U^A^U + U^d^U (1) 

whereas the physical observables are invariant under proper gauge transformations. The latter has to 
be everywhere smooth and it has to decrease at infinity in a suitable way (see, for instance, [I]). 

Since it is still unknown how to use in practice gauge invariant variable in Yang-Mills casqj, the 
most convenient practical choices for the gauge fixing are the Coulomb gauge and the Landau gauged). 
The standard approach to fix the gauge and to use Feynman expansion around the trivial vacuum 
An = provides one with excellent results when the coupling constant is small. 



In the cases of topological field theories in 2+1 dimensions [2] this goal has been partially achieved. 
2 The axial and light-cone gauge fixing choices are affected by quite non-trivial problems (see, for instance, [3]) 
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As it was for the first time discovered in [1] , a proper gauge fixing is not possible globally due to 
the presence of Gribov copies^. In the QCD case, this effect is very important in the non-perturbative 
regime. Even if one chooses a gauge free of Gribov copies, the effects of Gribov ambiguities in other 
gauges generate a non-perturbative breaking of the BRST symmetry |17j . 

It has been suggested to exclude from the domain of the path-integral gauge potentials which 
generate zero modes of the FP operator (see, in particular, [5] [7J [5] [S] [TU] [H]; two nice reviews are 
|12j |16|). In this framework, which is called the (refined) Gribov- Zwanziger approach to QCD, The 
region f2 around A^ = in which the FP operator is positive is called Gribov region: 

n d = f { Ap] d^A^ = and det &*D (A)^ > oj , (2) 
D(A)^ = d, + [A,,} , (3) 

where D {A)^ is the covariant derivative corresponding to the field A^. In the case in which the space- 
time metric is flat and the topology is trivial this approach coincides with usual perturbation theory 
when the gauge field A^ is close to the origin (with respect to a suitable functional norm 111]). At the 
same time, this framework takes into account the infra-red effects related to the partial elimination 
of the Gribov copies [8] [IB] |19j . When one takes into account thepresence of suitable condensates 
[2Q| [21] [22] [231 |25] the agreement with lattice data is excellent^ [2J] [25]. Within the (refined) 
Gribov- Zwanziger approach to QCD, the copy-free neighborhood of the trivial vacuum A^ = has to 
be identified with perturbative region of the theory. In this way, within the perturbative region, the 
standard recipe to fix the gauge and perform the perturbative expansion makes sense thanks to the 
lack of any overcounting. The appearance of copies of the vacuum satisfying the strong boundary 
conditions is one of the worst pathologies of perturbation theory. 

On the other hand, it is a well established fact by now that Yang-Mills theory may have knotted 
excitations [27] [28] [29] (the simplest non-trivial example corresponding to a closed tube: the "donut" 
or unknot with topology S 1 x D 2 ). Moreover, as it is well known (see, for instance, [36]), non-trivial 
topologies such as R x T 2 and MxT 3 are very important to understand the infinite volume limit 
of Yang-Mills theory as well as chiral symmetry breaking generated by the Casimir force (a related 
reference is [37]). In order to describe these situations one has to be able to define consistently 
Yang-Mills theory inside bounded regions with the non-trivial topologies. 

The main goal of the present paper is to analyze whether or not Gribov copies of the vacuum can 
actually exist inside a space-time region of the topology such as S 1 x D 2 (but more general topologies 

'furthermore, it has been shown by Singer [5], that if Gribov ambiguities occur in Coulomb gauge, they occur in all 
the gauge fixing conditions involving derivatives of the gauge field (See also [6]). Other gauge fixings (such as the axial 
gauge, the temporal gauge, and so on) free from gauge fixing ambiguities are possible but these choices have their own 
problems (see, for instance, [3])- 

4 The condition to have a positive Faddeev-Popov operator is not enough to completely eliminate Gribov copies in 
the Coulomb and Landau gauges, ft can be shown [TT] that there exist a smaller region (called the modular region) 
contained in the Gribov region which is free of gauge fixing ambiguities. However, it is still not clear how to implement 
the restriction to the modular region in practice. 

5 Recently, in [32], some criticism of this approach has been proposed. However, this argument is only valid in 
the perturbative framework (for instance, the gluons are assumed to be asymptotic states) while the changes of the 
Gribov- Zwanziger approach are non-perturbative in nature. 
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relevant, for instance, for the infinite volume limit of Yang- Mills as well as lattice QCD-such as T 3 and 
IR x T 2 , K x T 3 as well as the interior of an ellipse-will be considered). The absence of Gribov copies 
of the vacuum (satisfying strong boundary conditions) is a necessary condition for the existence of a 
perturbative region in the functional space around = 0. It will be shown that such requirement 
implies non-trivial restrictions on the possible sizes and shapes of the corresponding regions. In a sense, 
these results are more surprising than the ones obtained in [30J [31] [32J in which it has been shown 
that the pattern of appearance of Gribov copies strongly depends on the space-time metric. All the 
examples of non-trivial copies of the vacuum will be in flat metrics (but with non-trivial topologies). 
The present analysis will also show that the Gribov phenomenon strongly depends on the shapes and 
sizes of the bounded regions where one wants to study gauge theories. 

Unlike the equation for the zero-modes of the Faddeev-Popov operator, the issue of the appearance 
of Gribov copies of the vacuum is non-linear in nature and therefore can be quite complicated when 
non-standard topologies are considered. The technical tool necessary in order to analyze such an issue 
is a novel (to the best of authors knowledge) generalization of the hedgehog ansatz beyond spherical 
symmetry for the non-linear sigma-model. Such a generalization is quite interesting in itself since it 
allows to reduce the non-linear system of coupled partial differential equations corresponding to the 
equations of motion of the non-linear sigma model to a single non-linear partial differential equation 
which can be analyzed with the tools of solitons theory. 

The paper is organized as follows. In section two the notion of strong boundary conditions with 
non-trivial topologies will be analyzed. In section three, the relations between the Gribov copies 
equation and the non-linear sigma-model will be discussed. In section four, the generalized hedgehog 
ansatz will be constructed. In sections five and six various non-trivial examples of copies of the vacuum 
will be described. Some conclusions will be drawn in the last section. In the Appendix a novel way to 
implement both spherical and elliptical symmetries (which could be useful in the context of calorons) 
will be presented. 

2 Strong boundary conditions with non-trivial topologies 

Here, we will first review the definition of strong boundary conditions in the case in which the metric 
is flat and the topology trivial. In the present paper we will mainly analyze the SU(2) case but many 
of the present results also extend to other Lie groups. A useful starting point is the definition of 
non-Abelian charges Q a (see [1] [I] [38], [39]; see, for a detailed review, [10]): 



where M is the constant-time hypersurface in four-dimensional Minkowski space, is the two- 
dimensional sphere of radius R, Hi is the outer pointing unit normal to the indices i, j, k refers 
to space-like directions while a, b, c,...are internal SU{2) indices. 

Thus, one can define a proper gauge transformation U as a smooth gauge transformation which 




(4) 
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does not change the value of the charge as a surface integral at infinity: 




R 



d 2 Z (niE ia ) = lim f 
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) 



(5) 



Therefore, U is proper if it is smooth and approaches to an element of the center of the gauge group 
at spatial infinity. 

A Gribov copy U on a fiat and topologically trivial space-time satisfies the strong boundary condi- 
tions if U is proper. A copy of this type is particularly problematic since it would represent a failure of 
the whole gauge fixing procedure. Indeed, if the vacuum = possesses a copy fulfilling the strong 
boundary conditions, not even usual perturbation theory leading to the standard Feynman rules in 
the Landau or Coulomb gauge would be well defined. 

2.1 Strong boundary conditions on S 1 x D 2 

The simplest topology corresponding to a closed knotted tube is T = S 1 x D 2 : this case is relevant in 
the analysis of non-trivial topological excitations of glue-balls. In particular, there is a sound evidence 
supporting the existence of excited glue-ball states with the topology of T = S 1 x D 2 (see, in particular, 
[2"B] [25]). In this case, the spatial metric describing the region T reads: 



where the coordinate 4> corresponds to the S circle of T, while the coordinates r and 9 describe the 
disk D 2 of T. 

The radius of the disk is R and the boundary of T is: 



where Z 2 is the centeqj of SU(2). In other words, a smooth gauge transformation U = U(4>,r,9) is 
proper if 



ds 2 = [d(j)f + (dr 2 + r 2 d9 2 ) , 



(6) 



0<</><2vr, 0<6»<2vr, 0<r<.R 



37 = S 1 x 3D 2 . 



Thus, in this case a smooth gauge transformation U is proper if 



U(<t> + 2mvr, r, 9 + 2nvr) = U(<p, r, 9) , : 
U{4>,R,9) eZ 2 ,V<j>, 9 . 



m,n € Z 



(7) 
(8) 



Obviously, the center of 5(7(2) is made of two elements: ±1. 
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2.2 strong boundary conditions on T 3 

In the cases in which the spatial topology is T 3 (which, for instance, is relevant in the case of lattice 
QCD), the flat spatial metric describing T 3 reads: 

3 

ds 2 = >H (d&f , Xi € R , (9) 

8=1 

< <j>i < 2vr , 

where the coordinate (pi corresponds to the i-th factor S 1 in T 3 while Aj represents the size of the 
i-th factor S 1 in T 3 . In this case, due to the fact that <9T 3 = 0, a smooth gauge transformation 
U = U{<j>\ , 4>2 , 4>z) is proper if 

U((f)i + 2mi7r, </> 2 + 2m 2 7r,03 + 2m 3 7r) = U ((pi, <j>2, <P?i) , rrii e Z . (10) 

2.3 strong boundary conditions on IR x T 2 

The cases in which the spatial topology is R x T 2 may correspond to situations in which the bounded 
region in which we want to define perturbative Yang- Mills theory is much longer in one direction (the 
x axis) with respect to the other two. 

It is worth emphasizing here that, in the Euclidean theory, the topology R x T 3 is also very 
important in order to understand both the infinite volume limit in QCD (see, for instance, [36 ) 
and effects of chiral symmetry breaking in QCD (see, for instance, |37J). In particular, it has been 
shown in |37] the relevance of the Casimii0 force to understand chiral symmetry breaking. Hence, the 
importance of finite-size effects to explain these phenomena and the necessity to have a well-defined 
perturbation theory in such bounded regions with non-trivial topology require a deeper understanding 
of the Gribov problem in bounded region as well. As it will be shown in the following, the present 
results on the Coulomb gauge can be trivially extended to the case of the Landau gauge defined on 
R x T 3 . Thus, the gauge-fixing pathologies arising in this case have to be carefully taken into account. 

The flat spatial metric describing R x T 2 reads: 

2 

ds 2 = dx 2 + Y^ Xj (d(f>i) 2 , Xi E R , (11) 
i=i 

< < 2tt , 

where the coordinate (pi corresponds to the z-th factor S" 1 in T 2 while Xi represents the size of the i-th 
factor S 1 in T 2 . In this case, 9(lx T 2 ) is non-trivial and can be identified with the two limit point 
x — > ±oo of R. Thus, a smooth gauge transformation U = U(<pi, 4>2,x) is proper if 

7 Indeed, Casimir force is a genuine finite-size effect which, obviously, is not visible if the theory is analyzed in 
unbounded regions. 
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U{(f>i +2mi7r, 02 + 2m27r,x) = U(4>x,(/)2,x) , rrii G Z , (12) 
lim U((f>i,(f)2,x) £ Z 2 , V (pi, 4>2 ■ (13) 

3 Gribov copies and non-linear sigma model 

In the present paper we will analyze the copies of the vacuum in the Coulomb gauge but many of the 
present results can be easily extended to the Landau gauge. Since the Gribov copies of the vacuum in 
the Coulomb gauge can be seen as the Euler-Lagrange equations corresponding to a non-linear sigma 
model, we will first review some basic features of this model. The nondinear sigma model Lagrangian 
in D spacedike dimensions can be written in terms of group valued scalar field U. In the following we 
will consider the SU(2) group, so that the corresponding Lagrangian is: 

S = T / ^d D xTr [R l Ri] , (14) 
Ri = U^diU , (15) 



Ri — -Rj t a , 

t a t b = _ 5 ab 1 _ £ abc t c ^ 

[Ri,Rj] = CabR-iR'j , C c ah = —2e ca b , 

-.abc-.mnc ( cm rn ^X m \ 

£ £ — \°a °b ~ °a°b I J 

where D is the number of space dimensions, k is the sigma-model coupling constants, g is the determi- 
nant of the spacedike metric, 1 is the identity 2x2 matrix and t a the generator of SU (2) (where the 
Latin letters a, b, c and so on corresponds to group indices). The equation characterizing the appear- 
ance of Gribov copies of the vacuum in the Coulomb gauge (which corresponds to the Euler-Lagrange 
equations of the action in Eq. (HH)) read: 

= V* (U^diU) = , (16) 

where V 2 is the Levi-Civita covariant derivative corresponding to the metric (which, in the present 
case will have one of the forms in Eqs. ©, © and (fTT]) ). The case of the Landau gauge is analogous. 
The following standard parametrization of the SU(2)-vslued functions U{x l ) is useful: 

17(2:*) = yDl + Yat a , U~ 1 {x i ) = Y°l - Y a t a , (17) 
y o = Y °(x i ) , Y a = Y a {x i ) , (18) 
(Y°) 2 + Y a Y a = 1 , (19) 
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where, of course, the sum over repeated indices is understood also in the case of the group indices (in 
which case the indices are raised and lowered with the flat metric 5 ab ). Therefore, the Ri in Eq. (|15p 
can be written as follows: 

R c = £ abc Yad , Yb + Y°diY c - Y c diY° . (20) 



4 Generalized Hedgehog ansatz 

Due to the intrinsic non-linear nature of the Gribov copies of the vacuum, it is necessary to introduce 
a suitable technical tool which allows to study such phenomenon with non-trivial topologies. 

Here we will first discuss the geometrical interpretation of hedgehog ansatz as an effective tool to 
reduce the field equations of the non-linear sigma model (which, generically, are a system of coupled 
non-linear partial differential equations) to a single scalar non-linear partial differential equation. 
The geometrical analysis of such very important feature of the hedgehog ansatz allows to construct 
the natural generalization of the hedgehog ansatz for the non-linear sigma model. This section is 
interesting in itself since, to the best of authors knowledge, this is the first systematic reduction of the 
equations of motion of the non-linear sigma model (which is a system of coupled non-linear partial 
differential equations) to a single non-linear scalar equation beyond spherical symmetry. In terms of 
the group element U the standard spherically symmetric hedgehog ansatz reads 

U = 1 cos a(r) + n a t a sin a(r) , U" 1 = 1 cos a(r) — n a t a sin a(r) , (21) 
n 1 = sin 9 cos <fi , n 2 = sin 9 sin <f> , n 3 = cos 9 , (22) 

where r, 9 and (f> are spherical coordinates of the flat Euclidean metric, the n a are normalized with 
respect to the internal metric S a f,. In terms of the variables 1"° and Y a the hedgehog ansatz corresponds 
to the following choice: 

Y° = cos a(r) , Y a = n a sin a(r) , 5 ab n a n b = 1 . (23) 
Thus, the expression for Rf in Eq. (|20p reads 

RI = 5 dc (sin 2 a) e abd n a din b + n c d ia + sm ^ g-H c . (24) 

The equations of motion for the non-linear sigma model in Eq. (]16p corresponding to this ansatz read 

(V*fifca) +£sin2a = , (25) 

where L is a suitable function of the radial coordinate r (see below). Indeed, this is a very important 
and non-trivial characteristic of the spherically symmetric hedgehog ansatz which reduces a system 
of coupled non-linear partial differential equations to a single scalar equation. It is easy to see that 
in the case of the flat three-dimensional Euclidean metric on M 3 in spherical coordinates the above 
reduces to the standard Gribov pendulum of the vacuum (see [1] [T2]). 
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The previous analysis suggests to analyze which are the geometrical conditions which allow to 
reduce a system of non-linear coupled Partial Differential Equations to a single scalar PDE. Let us 
consider the following generalization of the hedgehog ansatz: 

U = 1 cos a + n a t a sin a , U^ 1 = 1 cos a — n a t a sin a , (26) 

5 ab n a n b = 1 , a = a (x*) , n a = n a (jc*) , (27) 

such that 

(did) (V*n a ) = , (28) 

(V*Si) n c = 2Ln c , (29) 

where L (which has to be the same for all the non-vanishingf§ n c ) may depend on the space-time 
coordinates. 

The important point is that, when the conditions (|28p and (|29p are satisfied, the function a (which 
in the usual spherically symmetric hedgehog ansatz can depend only on the radial coordinate r) can 
now depend on any set of coordinates and the functions n a (which in the usual case have to coincide 
with the unit radial vector) can be adapted to cases in which there is no spherical symmetry (and, 
therefore, no natural radial coordinate). Indeed, the expression for Rf in Eq. (|20p reads 

R\ = 5 dc (sin 2 a) e abd n a din b + n c ^a + Sm Ofi . (30) 

Hence, provided Eqs. (128p and (I29p are fulfilled and taking into account that the expression 

{V l n a ) (d % n b ) = gV (d^) [d,n b ) 

is symmetric under the exchange of a and b (so that its contraction with e a M vanishes) the coupled 
system of non-linear Partial Differential Equations (PDE) corresponding to the equations of motion 
of the non-linear sigma model reduces to the single scalar non-linear partial differential equation of 
sine-Gordon type: 

(V^a) +L sin 2a = . (31) 

Indeed, it is easy to see that in the usual spherically symmetric flat ansatz in Eqs. (|2ip and (|22D 
satisfies the conditions in Eqs. (I28p and (j29[) with L = —1/r 2 and that, correspondingly, Eq. (|31|) 
reduces to the usual Gribov pendulum (see [1] [E]). However, it is worth noting here that the present 
derivation is far simpler (just three lines of computation) and, furthermore, it discloses in a very clear 
way the geometry behind the hedgehog ansatz. 

The usual Gribov pendulum equation on flat and topologically trivial three-dimensional spaces 
(see [1] [12]) does not coincide with the sine-Gordon equation due to the radial dependence of the 
determinant of the three-dimensional Euclidean metric in spherical coordinates. For these reasons, 

8 Note that, in order for the solution to be non-trivial, at least two of the three n c have to be non- vanishing. 
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there are no analytic examples of copies of the vacuum in this case. On the other hand, using 
the present generalized hedgehog ansatz, one can construct many analytic examples of copies of the 
vacuum. These examples are very useful since they allow to express the requirement of absence of 
copies of the vacuum in terms of explicit constraints on the shapes and sizes of the corresponding 
bounded regions where one wants to define a perturbative region. 

4.1 Regularity conditions 

Another important requirement to satisfy the strong boundary conditions is the regularity of the pure 
gauge field U~ 1 dU everywhere. In the standard spherically-symmetric case, the only problematic 
point is the origin of the coordinates in which the spherical coordinates system is not well defined. 
One way to understand this condition in general is to analyze the behavior of the one-form 

R c = Rl A dx i (32) 

(where R\ is given by (|30p and a satisfies (|3ip ) when one changes coordinates from the Cartesian 
to a coordinates system which is singular somewhere. Indeed, in the following we will need to use 
non-Cartesian coordinates systems which are adapted to the symmetry of the problem (such as the 
elliptic coordinates systems in the next sections). 

In the case of Cartesian coordinates, the regularity at the origin can be read by looking at the 
components of the one-form i?? since the coordinates system is regular everywhere. On the other 
hand, if one is in a spherical or elliptical coordinates system, the angular coordinates are not well 
defined at the origin and some extra care is required. For instance, to disclose the singularity of the 
one forms d9 and d<f> at the origin one can analyze the Jacobian of the transformation from Cartesian 
coordinates (x n ) to spherical coordinates (x 3 ): 

dx H = J)dx j . (33) 

Due to the regularity of the Cartesian coordinates, the Jacobian itself encodes the informations of the 
singularity at the origin. In the spherically symmetric case, for example, if one takes into account 
the Jacobian, the regularity condition at the origin coincideqj with the usual one derived in Cartesian 
coordinates [12]. 

The requirement for regularity can be described as follows. Let's call x s a point in which the 
coordinates are singular and let's consider the following behavior for a when x — > x s 

a(x) nvr + /3/(x) (34) 

X— }X a 

9 In the spherical case the factor sin# which appears in the denominator of many components of the Jacobian is canceled 
out by a similar factor in the numerator of the spherically symmetric hedgehog ansatz in the normalized internal vector 
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with / (x) — > and f3 an arbitrary constant. Then, a copy of the form (|30p satisfy 

RC ^ /3 2 § dc f 2 EaUn a d . n b + ^ n c d j + ^ fd . n c (35) 

Then using (j33[) we see that f(x) must be such that 

R\ (J" 1 ); (36) 

(where (J -1 )* is the inverse Jacobian of the transformation from Cartesian to the curvilinear coordi- 
nates system of interest) are regular functions in x = x s . 

5 Explicit examples of vacuum copies in non-trivial topologies 

Here we will discuss many example of vacuum copies. 
5.1 S 1 x D 2 topology 

Here we will consider the following ansatz corresponding to the metric in Eq. (J6j) : 

a = a(r) , n 1 = cos (9) , n 2 = sin (9) , n 3 = . 

It is easy to see that the above ansatz satisfies the conditions in Eqs. ([28]) and ([29]) with the following 
L: 



n 1 = -^n 1 , (V*$) n 2 = -^n 2 

Thus, in this case the equation ()31|) for the copy of the vacuum reduces to 

rd r (rd r a) — 7 sin 2a = 44> (37) 



d 2 a 1 

7 sin 2a , 7 = - , r — tq = log r , (38) 



Q T 2 , , 2 



E= l - 
2 



fda\ 

I — I + 7 cos 2a 



(39) 



r-r = ± y (40) 

Jq(to) V 2 ^-7cos2y 

where ro and are integration constants. 

In terms of the radial coordinate r one has to require that the copy is regular at the origin and 
that it approaches an element of the center when r = R which is the boundary of D 2 . Without loss 
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of generality one can require (see Eq. ([26)) ): 



a — > , a(R) = mr , n £ Z . 

r— >-0 



(41) 



In terms of the variable r the above conditions read 



a — > , a — > mr 

(t— to)— >■— oc (r— to)— »t* 

r* = log i? . 



(42) 



Hence, we have to analyze under which conditions on the parameter 7 and on the integration constants 
E and tq Eq. (f39|) can have solutions satisfying the conditions in Eq. (f4U|) . The simplest way to answer 
to this question is to interpret Eq. (|39j) as the energy conservation of the following one-dimensional 
problem in which r plays the role of the effective time : 



V(A) = 27 cos A , A = 2a , E = 4E . 



(43) 
(44) 



The boundary conditions in Eq. (|42p in terms of ^4 read (we will consider the case n = 1) 



^ , A ->• 2vr 

(t— to)— >— 00 (t— to)— >t* 



(45) 



The above conditions mean that and 27r have to be two turning points and that the corresponding 
period has to diverge: 

dy 



1 



2~ 



cos y 

In order to satisfy this condition, it is enough to choose 



00 



(46) 



E = 27 = 1 



(47) 



the explicit form for the copy being 



a = 2 arctan 



r 
fo 



(48) 



with fg an integration constant (whose role will be described in a moment). 

In terms of the original radial variable r, the condition in Eq. (|4ip can only be satisfied"^*! if 
R — > 00. In physical terms, this means that the condition in Eq. (|4ip can be fulfilled when R is very 
large compared to fo: 

R > f . 

Obviously, in the present case the only two natural lengths are the radius of the disk R and the 



10 The reason is that the copy profile in Eq. (|48|) is an increasing function from to 00 and its maximum is ty. Note 
that, however, the copy profile approaches the asymptotic value very rapidly. 
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perimeter of the S 1 factor (which has been set to 2tt) of the donut S 1 x D 2 and consequently fg 
represents the size of S 1 . Thus, the previous analysis tells that when the radius of the disk is much 
larger than the perimeter of S 1 the donut is on the verge of supporting smooth copies of the vacuum 
satisfying the strong boundary conditions. Hence, in order to avoid this pathology the donut cannot 
be too "fat". 

It is easy to check that such copy is regular at the origin. The components of U~ 1 dU in Cartesian 
coordinates for this case are given by 

Ri{J~% , i = r,9 , j = 1,2 

With 



J- 1 



For r — > copies must decay as 



cos(#) sin(#) 

sin(fl) cos(fl) 



a -»• /3r , (49) 

r->0 



for some f3 S R and the profile in Eq. (|48|) satisfies this condition. 
5.2 T 3 topology 

As it has been already emphasized, the T 3 topology is relevant, for instance, in the analysis of the 
thermodynamical limit. Here we will consider the following ansatz corresponding to the metric in Eq. 
©: 

a = a(4>i) , n 1 = cos (pfo + #3) , n 2 = sin (pfa + qfo) , (50) 
n 3 = , p,q£Z . 

It is easy to see that the above ansatz satisfies the conditions in Eqs. (|28p and (|29p with the following 



Thus, in this case the equation (|3Tj) for the copy of the vacuum reduces to a flat elliptic sine-Gordon 
equation: 

d 2 \ 

—2 a = 7 sin 2a , (51) 



2 \ \\2j V A 3 
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In order for a solution of Eq. (|5ip to define a copy of the vacuum satisfying the strong boundary 
condition (see Eqs. (f2Uj) and ([TO]) ) it is necessary to require 



a(4>i + 2to7t) = a(0i) , m £ Z , (53) 



where we will consider m = 1 in the following. As in the previous subsection, Eq. (|51|) can be reduced 
to a first order conservation law: 



E = l 
2 



± 



/ da x 2 



+ 7 cos 2a 



l, 

d y 

a(4> ) V 2E - 7 cos 2y 



(54) 



where 0o and £7 are integration constants. Hence, we have to analyze under which conditions on the 
parameter 7 and on the integration constants E and (j)Q Eq. (|54p can have solutions satisfying the 
conditions in Eq. (|53|) . As in the previous subsection, it is possible to interpret Eq. (|54|) as the energy 
conservation of a one-dimensional problem in Eqs. (|43p and (|44p in which 0i plays the role of the 
effective time t = (pi. However, in this case the boundary conditions in Eq. (|53|) become: 

A(t + 2tt) = A{t) . (55) 

Thus, given two consecutive turning points Aq and Ai\ 

E 

Aq = arccos — , A\ = 2ir — Aq , 
2 7 

one has to require that the the time needed to go from Aq to A\ is half of the period in Eq. (|55p . 
namely: 

dy 



T - T = / - _ = 7T . 

lA o J2(E - 27 cosy) 

The time for a particle to go from Aq to A\ runs from to infinity as E runs from —27 to 27. This 
means that, given a 7, there is always an E such that r — tq = n. Then, for —27 < E < 27, it is 
always possible to construct copies of the vacuum satisfying the strong boundary conditions. 
Furthermore, it is easy to see that the norm of the copy is finite: 



\U- 1 dU\\= j d 3 x^jtr (U~ l dU) 

JT3 



2 



^^A 1 ^/ Ta *((^) 2 + ^ a ((^) 2 + (^) 2 ))<». 

On the other hand, it worth emphasizing that while there is a common agreement on the importance 
of the strong boundary conditions, it is not clear yet whether or not it is mandatory to require the 
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finite norm condition. For instance, it is possible to construct configurations of the Yang-Mills gauge 
potential which have infinite norm but finite energy and/or action (see, for a recent discussion, [41] ). 

It is worth emphasizing that, at a first glance, one could think that the presence of Gribov copies 
of the vacuum prevents one from using the Gribov semi-classical approach whose aim is, of course, to 
eliminate Gribov copies from a suitable neighborhood of the vacuum itself. However, the copies which 
can be eliminated using the Gribov semiclassical approach are "small" copies, namely zero- modes of 
the Faddeev- Popov operator (which is a linear condition). On the other hand, the vacuum copies 
which have been constructed here are solutions of the full non-linear equation for the copies and they 
would disappear if one would only consider the linearized equation (which reduces just to the Laplace 
equation in all the cases discussed in the present paper). Therefore, the Gribov semiclassical approach 
can be applied, for instance, to the R x T 3 or T 4 cases. For instance, in the T case in the Landau 
gauge, the inverse Gribov propagator (A^)" 1 would have, schematically, the usual form 

(Acr 1 = U + 1 U~ 1 , 

where 7 is the non-perturbative Gribov mass parameter determined by the usual gap equation [3] and 
□ is the Laplacian on T 4 . Hence, the only technical difference with respect the usual case would be 
that, in order to invert the Laplacian one should use the Fourier series instead of the Fourier transform. 
Indeed, in this way one can eliminate, from a suitable neighborhood of the vacuum = 0, zero modes 
of the Faddeev-Popov operator but the copies of the vacuum discussed above cannot be eliminated 
with this procedure. The physical consequences of this fact are very interesting also from the point 
of view of lattice QCD and are actually under investigation (for the relations between the continuous 
analysis of Gribov copies and the lattice see also [T3] , [J2] and [TH] ) . 

5.3 IxT 2 topology 

As it has been already explained, this topology is very important, for instance, in relation with the 
infinite volume limit of Yang-Mills theory as well as chiral symmetry breaking. 

Actually, as far as the infinite volume limit and chiral symmetry breaking are concerned (see, for 
instance, [36] [37]), the setting corresponds to the Euclidean theory defined on R x T 3 . Thus, in this 
case the Landau gauge is to be preferred. However, the present construction trivially extends to this 
situation as well. One has just to interpret the coordinate x in the metric in Eq. (jlip as the Euclidean 
time in order to apply the present analysis to this case. 

Here we will consider the following ansatz corresponding to the metric in Eq. (jlip : 

a = a(x) , n 1 = cos (p<f>i + q<f>2) , n 2 = sin (p0i + q<j) 2 ) , (56) 
n 3 = , p,q eZ . 

It is easy to see that the above ansatz satisfies the conditions in Eqs. (|28p and (|29p with the following 
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L: 



Thus, in this case the equation ([31]) for the copy of the vacuum reduces to a flat elliptic sine-Gordon 
equation: 



(—)<* 

\dx 2 J 



7 sin 2a , (57) 



7=^nf ) +m i . 



2 \ V X2J 



In order for a solution of Eq. (|57p to define a copy of the vacuum satisfying the strong boundary 
condition (see Eqs. ([26]) . (fl"2j) and (fT3|) ) it is necessary to require 

lim a(x) = rnir , lim a{x) = nir , m,n £ Z . (59) 

a;— >+oo x— 00 

Similarly to the previous subsections, Eq. (f5T|) can be reduced to: 



E = I 
2 



X — 



9a \ 

— — + 7 cos 2a 
ox J 

a(x ) V 2 ^ - 7 cos 2y ' 



(60) 



where and -E are integration constants. Hence, we have to analyze under which conditions on the 
parameter 7 and on the integration constants E and xq Eq. (|60p can have solutions satisfying the 
conditions in Eq. (|59|) for some m and n. As in the previous subsections, (|6U|) has the form of energy 
conservation of a one-dimensional problem (|43p with the identifications (|44p . In this case the strong 
boundary conditions in Eq. (|59p become: 

lim A(x) = 2rmr , lim ^4(x) = 2mr , m, n £ Z . (61) 

x— >+oo X— ^— 00 

For m = 1 and n = this means that and 2-7T have to be two turning points with infinite period. 
As in (|46p . this condition is ensured for E = 2-y. Furthermore, in this case as well it is easy to see 
that the norm of the copy is finite: 



- l dU\\= [ ^d 3 xtr(U- 1 dU) 2 < 

JT 2 xR 



\\l 

where 

a (x) = 2 arctan [exp (x — xq)] 
(xq being an integration constant) is the solution of Eq. (|60p . 
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6 Elliptic Gribov Pendulum 

In this section, we will analyze the elliptic generalization of the Gribov pendulum equation. It is easy 
to see that if one would try a naive generalization of the spherical Gribov pendulum [3] to elliptic 
coordinates (which reduces to the usual case in the spherical limit) then the system of equations for the 
Gribov copies of the vacuum do not reduce to a single scalar equation as it happens in the spherical case 
and, consequently, it would be extremely difficult to analyze the corresponding system of equations 
and the corresponding boundary conditions. On the other hand, within the present framework, one 
is lead to the two ansatzs in Eqs. ([6T?j) and ([76]) which do reduce the equations for the Gribov copies 
of the vacuum to a single scalar equation: this allows to study the strong boundary conditions in the 
usual way. To the best of authors knowledge, this is the first non-trivial elliptic generalization of the 
Gribov pendulum. 

This analysis is particularly important since it sheds light on how sensible the Gribov phenomenon 
is with respect to deformation of the spherical symmetry. Even if the usual Gribov pendulum equa- 
tio is analyzed in the unbounded region r G [0, oo[ [1] |12j . one implicitly assumes that similar 
results hold in a spherical bounded region. The reason is that, of course, experimentally gluons can 
live only within baryons and glueballs which are bounded regions. Therefore, it makes sense to study 
the arising of gauge copies of the vacuum in bounded region as well. Indeed, the results that no copy 
of the vacuum appear in the case of the usual Gribov pendulum also holds in the case of a bounded 
spherical region: as it will be now discussed, the elliptic case is more complicated. 

6.1 Prolate Spheroid 

The line element for of flat three-dimensional Euclidean space in prolate spheroidal coordinates is 
given by 



For a prolate spheroidal bounded region the coordinate ranges are given by < fi < R, v G [0, ir] , 
4> £ [0, 2it). The coordinate [i represents the elliptic radius since fi = const surfaces are ellipses. Such 
fj, = const ellipses have different eccentricities: large eccentricities (namely, large deformations from 
spherical symmetries) correspond to small \x while small eccentricities (namely, small deviations from 
spherical symmetry) correspond to large [i. Thus, in the above coordinates system, if one wants to 
analyze the limit of " large deformations from spherical symmetry" then one has to consider the small 
fj, region. On the other hand, if one wants to consider the almost spherical case, then the large \x limit 
has to be considered. As it will be shown in the following, the large deformation limit is the most 
interesting case. 

The Laplacian in this coordinates is given by 




(62) 



V 



l di 



1 




1 




a 



2 sinh 2 fi sin 2 v 



ii 



Actually, the whole issue of gauge copies is usually analyzed only in unbounded region. 
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Then, the following ansatz satisfies the conditions in Eqs. (|28p and (|29p 



a = a (fi,u) , n 1 = cos cj) , n 2 = sin , n 3 = , (63) 

with the following L: 



L 



2a 2 sinh 2 [i sin 2 z/ 

6.1.1 The strong boundary conditions 

Due to th fact that \x represents the elliptic radius, the strong boundary conditions in an unbounded 
region corresponds to the following conditions on a: 

a(n,v) -> nvr + f(fJ,,u) , (64) 

/x-S>0 

a(/j>,v) — > mir , n, m G Z . (65) 

where /(/x, ^) ensures regularity of the copies. On the other hand, if one is analyzing the theory in a 
bounded region of elliptic radius R then the condition in Eq. (|65p has to be replaced by 

a (//, u) — > mir , m £ Z , (66) 

since the boundary of the region is the surface fi = R and one has to require that the gauge copy 
belongs to the center of the gauge group on the boundary of the region itself. 

6.1.2 Regularity conditions 

In the prolate spheroidal case, the inverse of the Jacobian in the transformation (|33|) reads 



j-i = ( d(n,v,4>) \ _ 1 
d(x,y,z) 



/ cosh fj, sin v cos <f> cosh fj, sin v sin 4> sinh fi cos v \ 

5 5 — 5 5 — 5 5 — \ 

sinh ^t+sinh" !/ sinh /i+sinh f sinh ^t+sinh v 

sinh cos v cos sinh fj, cos 1/ sin </> cosh fj, sin 1/ 

5 5 — 5 5 — 5 5 — 

sinh ^t+sinh sinh /^+sinh 1/ sinh /i+sinh ^ 

sin cos <f> q 

\ sinh [i sin i> sinh sin v / 



(67) 



Thus, singularities could appear for fj, = and v = 0. Following the prescription in Eq. (f35j) . near 
points with \x = and/or ^ = copies behave as 

^ ->• /J 2 ^/ 2 ^^^ 6 + (3fd^n c 

where / must be such that functions (|36|) are regular. A sufficient condition which ensures regularity 
is: 
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Besides the singularity at ^ = 0, there are singularities for v = and u = ir as well (which are similar 
to the l/sin# singularity in the spherical case). However, in the elliptic case, the sinz^ factor is not 
automatically canceled by the internal vector n l of the hedgehog ansatz (which only depends on (j> in 
the present case). Therefore the profile function a(yu, v) has to take care of this divergence. 

6.1.3 Prolate pendulum 

Eq. ([3T|) in the prolate elliptic coordinates reduces to the following elliptic prolate Gribov pendulum: 

2 3 — [®u + ®u + co *h + c °t v ®v\ a = 2 5 — s ^ n ^ a ■ (^) 

sinh [X + sin v 2 sinh fj, sin v 

Now we consider the limits fx — > and /z — > oo which correspond to the cases of large and small 
deformation from spherical symmetry respectively. 



For n — > 0, equation (|68p takes the form: 

\d 2 Ll + d"l + coth fidn, + cot ^9^1 a = g — sin 2a , (69) 

2 sinh fi 

Interestingly enough, in this limit the ansatz a = a (fi) is consistent and the above equation 
reduces to 



da da 1 

— 77 + coth^— = „ 

dfi* an 2 sinh \jl 

With the following change of coordinate 



——77 + coth a— = k — sin 2a , (70) 



t = In tanh j (71) 

[i = 2 arcth e T (72) 

x 

sinh (arcth x) = —=== (73) 



the equation (l70|) can be written as 

d 2 a (t) 1 



dr 2 2 



sin 2a (r) . 



Of course, in this case, the natural boundary conditions correspond to the analysis within a 
bounded region in Eqs. (164)) and (|66l) in which, in suitable units, the radius R is very small 



#«1. 

In terms of the coordinates r in Eq. (|7ip Eqs. (|64p and (|66p 

a (r) — )• 2nir , 



-oc 
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a(r) 



2miT , n, m G 



In the limit of large prolate deformations, the equation can be integrated analytically: 



1 

E =2 



c£a \ 1 

"l - + - cos 2a 

dr I 2 



T - T 



ol{t) 



dy 



a ( T o) J2E- ±cos2y 



(74) 
(75) 



This result discloses in a very clear way how sensible the Gribov phenomenon is not only to 
topology but also to the shapes of the region where it is analyzed. Indeed, even if the present 
analytic solution cannot satisfy the regularity conditions in v = and v = ir (which are the 
north and south poles of the ellipse) since it is ^—independent, a small deformation of the ellipse 
at the poles could eliminate the necessity to require regularity of the solution at v = and v = tt 
and consequently could give rise to the sudden appearance of a copy of the vacuum. Hence, 
the present analysis strongly suggests that large prolate deformations from spherical symmetry 
can be quite pathological. These results can be relevant in the cases in which one is analyzing 
Yang-Mills theory in bounded regions. 



• For fi — y oo, equation ([68]) reduces to: 

[d 2 + <9 2 + <9 M + cot vdj\ a = - esc 2 v sin 2a . 

Since the large \i limit corresponds to small deviations from spherical symmetry, in this case one 
should recover the standard results on the absence of vacuum copies. 

6.2 Oblate Spheroid 

The line element for a flat three-dimensional Euclidean space in oblate spheroidal coordinates is given 
by 

ds 2 = a 2 (sinh 2 /j, + sin 2 v) {d[i 2 + du 2 ^ + a 2 cosh 2 fi cos 2 udef) 2 

For a oblate spheroidal bounded region the coordinate ranges are given by < \i < fi* , v € [0, tt] , (j) G 
[0, 2tt). Also in this case \x is the elliptic radius since (i = const surfaces are ellipses with eccentricities 
which decrease with [i. As in the prolate case, if one wants to analyze the limit of "large deformations" 
then one has to consider the small [i region. On the other hand, if one wants to consider the almost 
spherical case, then the large u limit has to be considered. 
The Laplacian in this coordinates reads 

V % di = 2 1 r -=- T [d 2 +dl+ tanh fjid^ + tan vd v \ + — — ^ —dl 

a A ( sinn fi + sin v) a z cosh [i cos z v 
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Thus, the following ansatz satisfies the conditions in Eqs. (|28|) and (|29p 



a = a (fj,, v) , n 1 = cos (j> , n 2 = sin 4> , n 3 = (76) 

with the following L: 



2a 2 cosh fi cos f 

6.2.1 Strong boundary conditions 

As in the prolate case, the strong boundary conditions in an unbounded region corresponds to the 
following conditions on a: 

a (//, v) — > 2mr , (77) 
a (/j,, u) — > 2m7r , re,m G Z . (78) 

On the other hand, if one is analyzing the theory in a bounded region of elliptic radius R then the 
condition in Eq. (fTSj) has to be replaced by 

a(/j,,u) — > 2m,7r , n,mGZ. (79) 

Regularity conditions As it has been explained in the previous subsection, the informations about 
the singularities at the origin are encoded in the Jacobian. Also in this case, besides the singularity 
for = 0, there are singularities for v = and v = tt as well Since the sin is factor is not automatically 
canceled by the internal vector n l of the hedgehog ansatz, the profile function a(fi, v) has to take care 
of this divergence: a sufficient condition to ensure regularity is: 

/(/i,i/) -> uV 

6.2.2 Oblate pendulum 

Hence, the equation for copies of the vacuum (|3ip in the oblate case now takes the form of the following 
oblate Gribov pendulum 

2 2 — i^u + ^ + tarA fid^ + tan vd u ] a = k — sin (2a) . (80) 

sinh /i + sin v 2 cosh fi cos 2 v 

Now we consider the limits fi — > and /i — > oo corresponding to large and small deformation from 
spherical symmetry respectively. 



For fi — > (f80|) reduces to 

[<9^ + <9 2 + tan vd^ a = — tan 2 sin (2a) 
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Unlike the prolate case, in the present case it is not possible to find analytic solutions in the 
limit of large deformations. 

• For fx — > oo ([SUj) reduces to 



Also in this case, since the large \i limit corresponds to small deviations from spherical symmetry, 
one should recover the standard results on the absence of vacuum copies. 

It is worth emphasizing that in the large deformation limit the prolate and oblate Gribov pendulum 
equations (|68p and (|80p differ significantly. In particular, unlike the spherical or oblate cases, in the 
prolate case the Gribov pendulum equation can be integrated exactly. This strongly suggests that 
prolate deformations from spherical symmetry are more pathological than oblate deformations. Indeed, 
the present results strongly suggest that the strong deformation limit of Eq. (|68p may support copies 
of the vacuum. This analysis appears to be quite relevant as far as the issue of gauge copies in bounded 
region is concerned. 

7 Conclusions and further comments 

In this paper the arising of Gribov copies in regions with non-trivial topologies (such as closed tubes 
S 1 x D 2 , or Ix T 3 ) but flat metric has been analyzed. The technical tool has been a generalization of 
the hedgehog ansatz beyond spherical symmetry. Such a generalization of the hedgehog ansatz is very 
interesting in itself since, in the case of the non-linear sigma-model, it provides one with a geometrical 
recipe to reduce the the field equations of the non-linear sigma model (which are a system of coupled 
non-linear partial differential equations) to a single scalar non-linear partial differential equation even 
when there is no spherical symmetry. This ansatz allows to construct many analytic examples of 
Gribov copies of the vacuum. Moreover, the elliptic Gribov pendulum has also been derived (to 
the best of authors knowledge, for the first time) both in the prolate and oblate cases. Our results 
suggest that large prolate deformations from spherical symmetry are likely to be more pathological 
than the oblate deformations. The requirement of absence of Gribov copies of the vacuum satisfying 
the strong boundary conditions implies geometrical constraints on the topology, on the shapes and 
on the sizes of the regions with non-trivial topologies (such as upper bounds on the deviations from 
spherical symmetry or constraint on the shape of the donut S 1 x D 2 ). Moreover, we have shown that 
in the case of a flat metric but with the topology of T 3 it is possible to construct copies of the vacuum 
satisfying the strong boundary conditions and with finite norm. The present results are interesting in 
relations with the infinite volume limit of Yang-Mills theory (which is related to the MxT 3 topology) . 
Indeed, one of the main point of the present analysis has been to show that when the T 3 (as well as 
T ) topology explodes to R? (or to R ) the vacuum copies disappear since neither. R s nor R 4 support 
copies of the vacuum. Hence, as far as the Gribov copies are concerned, such limit is not smooth and 




(81) 
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should be studied more carefully. This analysis is also relevant in all the cases in which gluons are 
confined in regions of finite sizes with non-trivial topologies such as in the cases of knotted flux tubes, 
lattice QCD and so on. Due to the close relation between Gribov ambiguity and confinement the issue 
of the Gribov copies in bounded regions (both on flat and on curved space-times) is very important 
and worth to be investigated. 
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Appendix: Landau gauge 

In this Appendix we will present two applications of the present generalized hedgehog ansatz which lead 
to a novel realization of spherical and elliptical symmetries respectively. The new way to implement 
spherical and elliptical symmetry corresponds to a configuration in which the internal vectors n l of 
the generalized hedgehog ansatz depend only on the Euclidean time. This realization of the spherical 
symmetry is only possible within the present "generalized hedgehog" framework and it could be useful 
in the context of calorons which are instantons which are periodic in Euclidean time. In both cases 
we will construct the Landau gauge pendulum. 

A novel spherical Case 

Let us consider the following line element corresponding to a four-dimensional euclidean space 

ds 2 = dr 2 + dr 2 + r 2 (d9 2 + sin 2 <9# 2 ) , < r < 2vr , 

where r plays the role of euclidean time (which is a periodic coordinate of period 2ir) and the spatial 
section is written in spherical coordinates. This situation is relevant in the cases in which one wants 
to describe finite temperature effects. An ansatz satisfying conditions ([2"8]) and (f29|) suited to deal this 
situation is 

a = a(r) , n 1 = cos (wr) , n 2 = sin (ur) , n 3 = 

with 
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Then, the equation for Gribov copies (|3ip takes the form 

2 L? 

d 2 a + -d r a = — sin (2a) , (82) 
r 2 

which obviously is not equivalent to the usual spherically symmetric Gribov pendulum [1] |12j . Defin- 
ing x = lnr we can write the above equation as 

to 2 

d 2 a + d x a = -y e 2x sin (2a) , 

another useful form is, defining y = —1/r, the following 

u? 

dyOt = ^4 sin (2a) 

The regularity condition at the origin can be analyzed as in the previous sections. 
Prolate Spheroidal Case 

Let us consider the line element of four-dimensional Euclidean space-time (in which the Euclidean 
time has been compactified to describe finite-temperature effects) in prolate elliptic coordinates 

ds 2 = (It 2 + a 2 (sinh 2 /x + sin 2 u) [dfi 2 + dv 2 ^ + a 2 sinh 2 \i sin 2 vd(f) 2 

An ansatz satisfying conditions (|28l) and ([29]) suited to deal this situation is 

a = a (/x, v) , n 1 = cos (ut) , n 2 = sin (wr) , n 3 = 



with 



L— 2l 



Then, Eq. (1311) reduces to 



1 r„9 „o , „ 1 W 2 



IS 2 + <9 2 + tanh/zeL + tan i/d„l a = — sin (2a) 
a 2 (sinlr fj, + sin 2 z/J ^ 2 

Also in this case, the regularity conditions can be studied as in the previous sections. 
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